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Recent  work  on  a  nonlinear  theory  for  gyro-devices  (i.e.,  gyrotrons,  gyro-twts, 
gyro-klystrons)  Is  described.  Differences  between  this  work  and  similar  work  (Thomas, 
Itftirkare  emphasized.  Comparisons  are  made  between  theoretical  predictions  and 
computer  simulations.  In  addition)  comparisons  are  made  between  theoretical  predictions 
and  60  GHz  gyrotron  experimental  data. 


A  NONLINEAR  THEORY  FOR  GYRO-DEVICES 


G.  B.  THOMAS 
VARIAN  ASSOCIATES,  INC. 

1.0  INTRODUCTION 

With  the  increasing  need  for  both  higher  power/higher  frequency  gyrotrons  and  for 
radar  quality  gyro-twts  and  gyro-klystrons,  it  is  important  that  a  more  complete 
understanding  of  the  basic  gyro-device  interaction  be  obtained.  It  is  no  longer  sufficient 
)ust  to  predict  power  and  efficiency.  Other  effects  such  as  harmonic  generation,  phase 
stability,  etc.,  must  be  well  understood  and  predictable.  To  accomplish  the  task  of 
modeling  such  nonlinear  effects,  a  nonlinear  theory  (Thomas  1981,  1983)  is  used.  Before 
work  can  begin  on  predicting  general  behavior  resulting  from  nonlinear  effects,  it  is 
important  that  predictions  of  the  nonlinear  theory  show  close  agreement  with  computer 
simulations  and  experiments.  The  nonlinear  theory  has  been  extended  to  the  point  where 
there  is  close  agreement  between  theoretical  predictions  and  work  (computer  simulation 
and  experimental)  done  at  Varian/Palo  Alto. 

In  this  paper,  the  recent  theoretical  work  based  on  the  concept  of  the  soliton  is 
described  in  detail  (SECTION  3.0).  Differences  between  the  previous  nonlinear  work  and 
this  recent  work  are  discussed.  These  differences  are  as  follows; 

1)  The  inclusion  of  the  effect  of  the  rf  fields  on  the  electron  equations  of 
motion. 

3)  The  inclusion  of  higher  cyclotron  harmonics  in  the  higher  order  dispersive 
terms. 

3)  A  more  complete  equation  for  the  number  density. 

4)  The  inclusion  of  dispersive  effects  in  the  direction  of  fast  wave 
propagation. 
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5)  The  inclusion  of  the  effects  of  the  rf  field  on  the  operating  range  of  the 
gyro-device. 

By  including  these  changes,  one  has  not  only  derived  a  nonlinear  partial  differential 
equation  which  describes  the  dynamic  evolution  of  the  nonlinear  interaction  between  the 
rf  fields  and  the  electron  plasma  in  a  gyro-device,  but  from  an  understanding  of  this 
equation,  one  can  determine  and  predict  the  saturation  mechanism  which  limits  the 
device's  operating  range.  This  saturation  mechanism  is  discussed  in  Section  3.0.  In 
Section  4.0  theoretical  predictions  of  the  magnetic  field  where  saturation  begins  and  the 
peak  rf  output  power  of  a  given  gyro-device  design  are  compared  to  available  simulation 
results  and  experimental  data.  Conclusions  and  areas  for  future  work  are  given  in  Section 
S.O. 

2.0  THB  NONLINBAR  THEORY 

The  general  technique  used  to  derive  the  nonlinear  partial  differential  equation 
describing  the  rf  field-electron  plasma  interaction  in  a  gyro-device  consist  of  the 
following  steps: 

1)  The  derivation  of  a  partial  differential  equation.  This  equation  is  derived 
using  the  method  of  characteristics  (Krall  and  Triveiplece,  1973,)  where  the 
effect  of  the  rf  electric  field  on  the  motion  of  the  electrons  is  included  in 
zero  order. 

2)  The  derivation  of  a  number  density  equation  describing  the  effect  of  the 
rf  electric  field  on  the  electron  beam.  This  equation  includes  the  effect  of 
the  dominant  nonlinearity  which  is  also  responsible  for  electron  bunching. 

3)  The  partial  differential  equation  and  the  number  density  equation  are 
coupled  together  and  result  in  a  nonlinear  partial  differential  equation. 
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4)  The  nonlinear  partial  differential  equation  Is  reduced  by  a  multiple  time 
scales  perturbation  technique  such  that  at  second  order  (i.e.i  slowly  varying) 
a  modified  nonlinear  Schrodinger  equation  (MNLSB)  is  obtained. 

5)  The  terms  of  the  MNLSB  are  ordered  such  that  If  the  coefficients  have 
the  proper  signs,  the  two  dimensional  cubic  nonlinear  Schrodinger  equation 
(NLSB)  is  obtained  at  zero  order. 

6)  The  effects  of  the  terms  that  modify  the  nonlinear  equation  from  the 
two  dimensional  NLSB  to  the  MNLSB  are  determined.  Typically,  additional 
dispersive  terms  result  in  growth  or  damping. 

7)  The  region  in  parameter  space  in  which  one  obtains  growing  soliton-llke 
solutions  is  determined. 

In  this  paper,  each  of  these  steps  is  decribed  in  sufficient  detail  to  make  the  general 
technique  clear.  When  possible,  reference  Is  made  to  previous  work,  therefore,  only  the 
most  recent  results  are  discussed. 

2.1  DERIVATION  OP  THB  PARTIAL  DIFFERENTIAL  EQUATION 

The  general  technique  for  deriving  the  partial  differential  equation  for 
gyro-devices  is  well  known  (CHU,  et.  al.,  1979).  The  Vlasov-Maxwell  equation  Is  solved 


for  a  perturbed  electron  distribution  function  using  the  method  of  characteristics.  For 
gyro-devices,  the  assumed  electron  distribution  function  is  typically  a  delta  function  in 
both  parallel  and  perpendicular  velocity  (i.e.,  fo*(6(P„)6(PA  -Pq)/2itPa  )).  In  previous 
work,  the  electron  equations  of  motion  only  included  the  effect  of  dc  fields.  However, 
wave  propagation  through  an  electron  plasma  can  be  significantly  different  if  the  wave 
drastically  alters  the  motion  of  the  individual  electrons.  Therefore,  the  effect  of  the  rf 
field  on  the  electrons  is  included  In  the  zero  order  equations  of  motion.  Including  the  rf 
field  in  the  electron  motion  is  the  first  major  difference  between  this  nonlinear  theory 
and  the  previous  work. 


Because  bunched  electrons  have  the  greatest  effect  on  rf  wave  propagation,  these 
are  the  electrons  of  interest.  The  assumption  Is  made  that  for  bunched  electrons  the  rf 
fields  appear  quasi- static,  it  should  be  noted  that  this  assumption  does  not  imply  that 
there  are  any  bunched  electrons  nor  does  this  assumption  indicate  anything  concerning  the 
mechanism  of  electron  bunching.  The  mechanism  of  electron  bunching  is  modeled  by  the 
number  density  equation  and  is  discussed  later.  The  validity  of  this  assumption  is 
supported  by  the  fact  that  the  angular  frequency  of  interest  is  no  longer  the  fast  time 
scale  parameter,  a,  but  instead,  is  the  slow  time  scale  parameter  t>-Q‘.  With  the 
geometry  shown  in  Figure  1  and  a  transformation  to  the  rotating  frame  of  reference,  the 
following  set  of  equations  is  obtained-, 


vx  =  -vx COS  (0-(o-Q‘)t)-Ey/Bo  (la) 

vy  =,  -v^  SIN  (Q-(u>-Q‘)X)  (lb) 

vz  =  v„  (lc) 

where  Q'-ft/y.  Notice  that  a  non- time  varying  BxB  drift  is  included  in  (la). 

Substituting  the  above  set  of  equations  into  the  Maxwell- Vlasov  equation  and 
integrating  over  velocity  space,  a  series  of  terms  are  obtained  which  contain  Bessel 
functions.  These  Bessel  funcions  are  treated  differently  in  the  present  nonlinear  theory 
than  in  both  previous  linear  and  nonlinear  theories.  In  linear  theory  (Chu,  et.  al.,  1979)  the 
Bessel  functions  are  used  to  describe  the  rf  field  character  of  the  various  modes.  In  the 
previous  nonlinear  theory  (Thomas,  1981),  expansions  of  the  Bessel  functions  for  small 
arguments  were  taken  to  obtain  higher  order  dispersive  terms.  These  expansions  assumed 
that  the  order  of  the  Bessel  functions  remained  at  one,  (l.e.,  N-l).  N  was  chosen  equal  to 
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one  because  the  assumption  was  made  that  the  gyro-device  interaction  Is  a  resonant 
Interaction  at  o-Q,  and  therefore,  any  harmonic  cyclotron  Interactions  are  negligible. 
This  assumption  is  not  valid  if  the  interaction  is  strongly  nonlinear  such  that  higher  order 
electron  motion  exists  as  a  result  of  the  interaction  of  the  rf  field  with  the  electron 
plasma.  Higher  order  electron  motion  is  no  longer  described  by  simple  motion  about  a 
magnetic  field  line,  but  rather,  contains  cycloiding  components  superimposed  on  the 
simple  motion  (Figure  2).  This  additional  motion  Indicates  that  the  electron  plasma  has 
become  weakly  turbulent  and  contributes  to  higher  order  dispersion.  Therefore,  Bessel 
function  expansions  with  N=»l,  N=2,  and  N=3  are  used: 


N2J^(a)  =  (l/4)(l2)a2-(l/16)(22la4*(l/64)(32)a6 


where  N=1  contributes  to  k°  terms,  N*2  contributes  to  k3  terms,  and  N=3  contributes  to 

aw  X 

4 

k^  terms.  The  inclusion  of  these  higher  order  cyclotron  harmonic  terms  in  the  higher 
order  dispersive  terms  is  the  second  major  difference  between  the  previous  nonlinear 
theory  and  the  work  presented  here. 

By  substituting  the  above  Bessel  function  expansions  into  the  various  component 
terms  resulting  from  the  integration  of  the  Vlasov-Maxwell  equation  over  velocity  space, 
the  following  equation  is  obtained; 


o  u  o  c  o  c.  O  E. 

A  +  B  +  c  --4-  +  D  -r4  -  0 
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(<*>  -  (k„v,(  +  n)/Y  -  kxEy/Bo)2 


(4c) 


_ (k„v„  +  n)/y _ 

(u  -  (k„v„  +  0)/y  -  kxEy/Bo) 


(4d) 


As  *  result  of  including  the  rf  electric  field  in  the  electron  equations  of  motion,  the  rf 
electric  field  appears  in  the  denominator  of  all  four  coefficients.  As  a  result  of  Including 
N*2(3),  terms  like  20(30)  appear  in  the  numerator  and  denominator  of  the  A(B) 
coefficient. 
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2.2 


The  number  density  equation  describes  the  effect  of  the  wave  on  the  electron 
plasma.  In  other  words,  the  number  density  equation  models  the  bunching  of  the  electrons 
which  is  a  result  of  the  ponderomotive  force.  This  ponderomotlve  force  is  not  the  usual 
ponderomotive  force  acting  along  the  direction  parallel  to  the  dc  magnetic  field,  but  is 
the  ponderomotive  force  acting  in  the  direction  of  the  rotating  electrons  resulting  in 
bunching  transverse  to  the  direction  of  wave  propagation.  The  previous  nonlinear  theory 
did  not  take  into  account  electron  motion  in  the  z  direction  when  considering  the  bunching 

A- 

mechanism.  Z  directed  motion  is  considered  in  the  more  complete  fluid  number  density 


equation  given  below; 


32n1  k  2  32n 

1  2  z 

3tz  Vz  kxz  ¥x2 


3x  ^  m  ' 


where  v  is  the  wave  group  velocity  and  n.  is  the  perturbed  number  density.  Assuming  a 
Z  1 

plane  wave  solution  and  transforming  into  the  rotating  frame  yields; 


k  2 

-(“  "  fi,)2  ni  +  V  if7  kx2  ni 

X 


JL  (  NL' 

3z  ^  m  ’ 


(6) 


By  substituting  the  expression  for  the  ponderomotive  force  into  (6),  the  following 
expression  is  obtained  for  n^; 
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2.3  DERIVATION  OP  THE  PULL  NONLINBAR  PABTIAL  DIPFBRBNTIAL  EQUATION 

The  full  nonlinear  partial  differential  equation  is  obtained  by  combining  the  two 

coupled  equations  (3)  and  (5)  into  a  single  nonlinear  equation.  (3)  is  coupled  to  (5)  because 

2 

it  contains  n.  in  its  coefficients  in  the  form  of  u  .  (5)  is  coupled  to  (3)  because  it  has 
1  P 

2 

|E|  in  the  ponderomotive  force  term.  Using  (7)  and  n*n  +n^,  the  following  is  obtained; 


+  r,  |e  |2e  =  o 
1  y  y 

2 

where  the  zero  superscript  indicates  n^n  in  a  as  follows; 

o  p 

A°  =■  A|  ,  etc.  (9) 
n»n„ 


and  where 


HL2  (1  -  D  °)  +  k  4  B°  -  k  2  C°~]  ) 

p  L—  X  Z  I  U) 


-  S2\2  ^z2  ' 


(10) 


The  above  technique  is  the  same  as  that  used  in  the  previous  nonlinear  theory. 


2.4 


Solving  the  full  nonlinear  partial  differential  equation  is  a  formidable  task.  Yet, 
using  a  perturbation  technique  based  on  expansions  in  small  amplitudes  around  an 
equilibrium  value  only  restricts  the  problem  requiring  one  to  first  linearize  (8). 
Linearizing  (8)  would  negate  all  the  effort  in  getting  (8)  in  the  first  place.  Therefore, 
instead  of  using  a  perturbation  technique  based  on  expansions  in  small  amplitudes,  a 
technique  based  on  expansions  in  time  scales  is  used.  This  is  the  multiple  scales 
perturbation  technique.  (Nayfeh,  1973). 

In  the  previous  nonlinear  theory  (Thomas,  1981),  the  following  ordering  of  partial 
derivatives  was  used; 


_3_ 

3t 


-iw 


6u 


c  3x' 


_3_ 

3x 


ik  + 
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3 

3x 1 


_3_ 

3y 


62 
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3y' 


_3_ 

3a 


6 


3 


3z ' 


(11a) 


(lib) 


(11c) 


(lid) 


where  6  is  a  small  parameter  indicating  the  order  of  time  variations.  In  other  words, 
higher  powers  of  6  indicate  more  slowly  varying  quantities.  It  should  be  noted  that  a 
traveling  wave  solution  is  assumed.  Variations  in  the  $  direction  (see  Figure  1)  are 
considered  second  order  since  the  y  direction  is  neither  the  direction  of  propagation  (x 
direction),  the  direction  of  bunching  (x  direction),  or  the  direction  of  velocity  drift  (the  z 
direction).  Notice  that  the  z  direction  is  considered  first  order  because  the  electrons 
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drift  in  this  direction.  Wive  propagation  along  the  z  axis  was  neglected.  Neglecting  wave 

propagation  along  the  I  axis  is  thought  to  be  a  weak  assumption;  therefore,  wave 

propagation  is  Included  in  the  analysis  by  including  Ik  in  (lid).  Assuming  that  ik  Isa 

z  z 

zero  order  (fast  time  scale)  effect,  (lid)  becomes*, 


al  *  lk2  +  5  jp"  (lid) 

Substituting  the  expressions  for  the  partial  derivatives  given  in  (lla)-(lld)  into  the 
nonlinear  partial  differential  equation,  (8),  one  obtains  the  following  set  of  equations: 


-A°k 


C°k 


-  D 


(12) 


t 


which  yields  the  phase  velocity  (v 


) 

pz 


* 


(13) 


i'  21„Uc  D°  ^  +  21V°  ^  -  41ltxV  £  +  2*0°  -  0 


(14) 


which  yields  the  group  velocity  (u  ), 

c 


2k  3B° 

(it  v 

uc‘k:  +  l^-  + 


fc°  ^z_  _  C°  3xH 

I  D°  ukx  Du  u  3z 1  I 


(15) 


tl  14  n°  3  1,  2„°  32  ^  „o  3  kz  32 

—  -21“D  2r'-2k’'“  3?7+c  [pr  -  0 
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(16) 


This  second  order  equation  models  a  slowly  varying  modulating  envelope.  The  cubic 

f 

nonlinear  term  is  also  considered  second  order  resulting  in  the  full  nonlinear  equation  as 
follows; 


(17) 

+  Tj  |Ey|2  Ey  *  0 

This  equation  is  a  modified  form  of  the  nonlinear  Schrodinger  equation  (MNLSE)  with  a 

cubic  nonlinearity.  It  is  clear  that  keeping  wave  propagation  in  the  z  direction  results  in 

the  additional  term  -(k  /k)(32E  /chr'3z').  If  k  =0,  the  above  nonlinear  partial  differential 

z  x  y  z 

reduces  to  the  equation  derived  in  previous  work.  This  additional  term  is  the  fourth  major 
difference  between  the  previous  nonlinear  theory  and  that  presented  here. 

2.5  ORDERING  THB  TERMS  OP  THE  MNLSE 

Unfortunately,  at  this  time  it  is  not  possible  to  solve  analytically 
the  MNLSE,  (17),  because  of  the  additional  spatial  variation  and  because  of  the 
dependence  of  the  coefficients  on  the  rf  electric  field.  Therefore,  it  is  necessary  to  treat 
the  additional  terms,  etc.,  as  perturbations  and  obtain  an  ordering  as  follows; 
zero  order  c° 

3E  ^2  r 

-2iwD°  (Ey  ■  const)  -r^  -  2kx2B°  (Ey  *  const)  (18) 

+  Tj  (Ey  -  const)  |Ey|2  Ey  *  0 


-2iuD°  _  2kx2B°  +  C° 


32Ey  32Ey 

3z'2  kx  3x'3z' 


I 
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By  setting  By  in  the  coefficients  equal  to  s  nonzero  constant,  the  time  varying  nature  of 
the  coefficients  is  eliminated  When  B°  is  negative  and  f  ^  is  positive,  (18)  becomes  the 
cubic  nonlinear  Schrodinger  equation, 


-i 


+ 


|Ey|2  Ey 


0 


(19) 


which  has  stable  (nonlinearly  stable)  soliton  solutions. 

*  i 

first  order  c 

C°  (Ey  *  const) 

t 

The  C°  term  is  considered  a  first  order  perturbation  to  the  zero  order  solution.  Higher 
order  effects  are  those  that  result  from  B  not  equal  to  a  constant,  but  where  B  in  the 

)  y  y 

coefficients  is  a  function  of  time. 

2.6  BPPBCTS  OF  THE  PERTURBATION  TERMS  ON  THE  NTLSB 

It  is  important  to  determine  the  effect  that  the  additional  term  (the  C°  term)  has 

► 

on  the  soliton  solutions.  Pollowing  Albowitz  and  Segur,  1979,  the  MNLSB  is  normalized  to 
the  following  form; 

IA  ♦XA-.+uA  +X|A|2A*0  (21) 

;  T  65  Tjrj 

where  uA.__  is  the  perturbing  term, 

un 

From  Albowitz  and  Segur,  1979,  the  following  results  are  obtained; 
a.  CASE  1  u>0.X>0 

* 

LA  +A-  -*A  ♦|A|2A»0  (22) 

*  9  5  tjti 


32Ev  k, 

-H?  -  5 


32Ey 

3x3z 


(20) 
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This  equation  is  unstable  to  symmetric  perturbations  (l.e.,  even  mode)  in  the  £  direction, 
b.  CASH  2  u>0,  \<0 


-‘VA{rVA,JA-° 


(23) 


This  equation  is  unstable  to  anti- symmetric  perturbations  (l.e.,  odd  mode)  in  the  £ 
direction. 

c.  CASB  3  u<  0.  \>  0 


IA  +A..-A  +|A|  A»0 

*  £4  Thl 


(24) 


This  equation  is  unstable  to  anti- symmetric  perturbations  in  the  £  direction, 
d.  CASB  4  u<0,  \<0 


-iA  +  A..+A  +|A|  A=0 

t  1 1  TPI 


(25) 


This  equation  is  unstable  to  symmetric  perturbations  in  £. 

To  summarize  the  above  cases,  the  additional  spatial  varying  term  results  in  the 
nonlinear  instability  of  the  soliton  solution.  The  type  of  perturbation  (symmetric  or 
anti-symmetric)  that  results  in  the  instability  depends  on  the  sign  of  the  additional  term. 
In  a  gyro-device,  the  nature  of  this  nonlinear  instability  is  related  to  the  bunching  of  the 
electrons  which  cause  the  perturbation  in  the  £  direction.  Depending  on  the  signs  of  the 
coefficients,  electron  bunching  around  the  zero  of  the  rf  field  (anti- symmetric)  or 
bunching  around  the  maximum  (symmetric)  of  the  rf  field  will  cause  the  nonlinear 
instability. 
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2.7  PARAMETER  BANGS  POP  GROWING  SOLITON-LIKB  SOLUTIONS 

The  final  step  in  the  general  technique  of  modeling  the  dynamic  nonlinear 
interaction  in  a  gyro-device  is  to  determine  the  region  in  parameter  space  where  one 
expects  growing  soliton-like  solutions  to  the  MNLSB.  In  order  to  determine  which  of  the 
possible  forms  (22)-(25)  of  the  equation  is  the  one  of  interest,  the  spatial  location  of  the 
bunched  electrons  must  be  considered.  Recall  from  section  2.2  that  the  ponderomotive 
force  acting  in  the  direction  of  rotation  is  responsible  for  electron  bunching.  The 
ponderomotive  force  causes  bunching  about  the  rf  electric  field  zero  which  is  an  anti¬ 
symmetric  perturbation.  Consequently,  only  cases  2  and  3,  (23)  and  (24)  remain.  The 
difference  between  these  two  cases  is  that  (23)  represents  a  system  where  the  envelope 
group  velocity  is  less  than  the  wave  phase  velocity  (see  (15))  and  (24)  represents  a  system 
where  the  envelope  group  velocity  is  greater  than  the  wave  phase  velocity.  In  the  case  of 
the  gyro-device,  the  velocity  of  the  bunched  electrons  is  less  than  the  wave  phase 
velocity  (in  the  x  direction,  not  in  the  z  direction),  so  (23)  applies.  To  obtain  (23)  from 
(17),  the  coefficients  must  have  the  following  signs; 

D°  Positive 
B°  Negative 
C°  Negative 
Tj  Positive 


The  range  of  d.c.  magnetic  field  values  over  which  growing  soliton  solutions  are 
obtained  is  determined  by  considering  the  magnetic  field  values  where  the  signs  of  the 
coefficients  change  from  those  listed  above.  First,  consider  the  B°  coefficient.  From 
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(4b),  B°  Is  negative  for  «<(k„v„+30)/Y-kx  ®y/B0*  The  of  interest  is  where  u  Is 

o  o  o 

near  Q  ,  so  the  B  is  always  negative.  Next,  consider  the  C  coefficient.  For  the  C 

term  (4c)  to  be  negative  the  following  is  obtained; 


(k,,v„  +  ~  k 


E  /B~~| 
x  y  oj 


2 

-oj  v 


1 


c‘7' 


(26) 


(26)  is  always  true  in  the  parameter  regime  of  interest.  Por  the  complete  C°  term  to  be 
negative 


32Ey  kz  32Ey 
3z~  y  kx  3x32 


(27) 


must  be  true.  From  (27),  the  following  constraint  is  obtained; 


(28) 


In  other  words,  for  large  values  of  k  ,  the  variation  of  the  field  in  the  3k  direction  must  be 

z 

small  or  the  variation  of  the  field  in  the  £  direction  must  be  large.  As  expected,  for 
positive  wave  numbers  (forward  wave  interactions),  the  slope  of  the  rf  field  is  such  that 
the  rf  field  strength  increases  along  the  £  direction.  Consider  the  D°  coefficient.  From 
(4d),  D°  is  positive  when 


Q  <  yw 


(29) 


with  B„ 
7 


assumed  equal  to  zero.  (29)  provides  the  constraint  on  the  maximum  allowable 
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magnetic  field  for  soliton  solutions.  The  [\  coefficient  provides  the  lower  magnetic  field 

constraint.  Notice  from  (10)  that  the  behavior  Is  dominated  by  the  1-D  term. 

l  o 

Therefore,  is  positive  when 


where  was  assumed  equal  to  zero.  Of  course,  for  the  case  where  =  0  the  simple 
analytic  solutions  given  in  (29)  and  (30)  are  not  complete.  However,  these  equations 
Illustrate  the  general  characteristics  of  the  operating  regime  of  the  gyro-device.  The 
complete  characteristics  of  the  operating  regime  are  given  in  the  next  section  which 
discusses  the  saturation  mechanism  of  the  gyro-device  interaction. 

3.0  THB  SATURATION  MECHANISM 

The  saturation  a  of  gyro-device  occurs  where  the  electrons  have  given  up 
sufficient  energy  so  that  soliton  solutions  to  the  nonlinear  partial  differential  equation  are 
no  longer  possible.  In  other  words,  the  maximum  magnetic  field  value  for  electrons  which 
have  lost  perpendicular  energy  to  the  wave  (B^(E^  *  0)|^  )  is  less  than  the  minimum 
magnetic  field  value  for  electrons  which  have  not  lost  any  perpendicular  energy  (Br  (E  » 

L  y 

0)|  ).  Therefore,  full  saturation  occurs  when, 

BL(Ey  =  0),v  “  BHfEy  “  0),v  /  (3l) 

The  constraint  in  (31)  allows  the  perpendicular  electron  energy  to  decrease  smoothly  from 
the  initial  value  to  the  final  value  without  interruption.  In  other  words,  growing  soliton 
solutions  are  possible  over  the  full  range  of  v^  values  given  a  proper  value  for  the  rf 
electric  field. 


r 


t 


f 


It  is  Importsnt  to  include  rf  field  effects  so  that  both  the  magnetic  field  for  peek 
output  power  (Bp)  end  the  ectuel  velue  of  the  output  power  et  Bp  cen  be  determined.  The 
rf  electric  field  (B^>  results  in  e  change  in  rotetionel  electron  velocity  through  en  Ey  X  Bq 
drift.  Cleerlyi  the  electrons  ere  eccelereted  or  decelereted  depending  on  the  sign  of  the 
rf  field.  The  sign  of  the  rf  field  depends  upon  the  locetion  of  the  electron  in  its  orbit. 

output  power  is  echieved  when  the  meximum  number  of  electrons  ere  loceted  in 
e  fevoreble  position  in  the  rf  field.  The  fevoreble  position  in  the  rf  field  is  defined  es  the 
spetiel  locetion  or  locetions  where  conditions  ere  such  thet  soliton  solutions  exist  to  the 
nonlineer  pertlel  dlfferentiel  equation.  Clearly,  meximum  output  power  would  be 
obtained  if  the  entire  electron  orbit  constitutes  e  fevoreble  position.  This  would  occur  if 
en  electron  in  either  the  positive  or  negative  rf  field  cen  contribute  to  the  growth  of  the 
field.  Therefore,  the  following  expression  is  obtained; 


» 


♦ 


VV 


.  ■  B„(B  -  >| 

Av  H  y 


v*« 


- 


(32) 


The  velue  of  magnetic  field  thet  satisfies  (32)  is  the  optimum  for  the  particular  design 

under  consideration.  The  Ay,  term  takes  Into  account  the  effect  of  the  E  .  X  li  drift  on 

rt  o 

the  rotetionel  velocity  of  the  electrons. 

Notice  thet  Av,  is  subtracted  from  both  v,  end  vA  ^  end  not  subtracted  from  v, 
end  added  to  vA  j  es  one  would  suspect  from  the  sign  of  shown  in  (32).  The  reason  thet 
-Av,  is  used  is  thet  v^  -  Ava  is  the  minimum  perpendicular  velocity  possible  from  the 
interaction  end  describes  electrons  which  have  lost  the  most  perpendicular  energy.  These 
electrons  ere  bunched  et  the  locetion  in  the  electron  orbit  where  the  rf  field  is  directed 
perpendicular  to  the  motion  of  the  electron,  (Figure  3). 

Having  obtained  the  optimum  Bp  end  the  meximum  perpendicular  energy  lost  to 
the  rf  wave,  the  rest  of  the  operating  characteristics  cen  be  determined.  First,  the 
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efficiency  Is  found  from  the  following. 


(yinitial  ~  1)  ~  (yfinal  ~  1) 
(yinihal  ”  1) 


(33) 


where  Is  used  to  obtain  end  va  {  -  is  used  to  obtain  YpiNAL  The  rf 

electric  field  strength  is  determined  from  Ava  as  follows! 


E,t 


(34) 


The  rf  electric  field  strength  is  converted  to  rf  output  power  by  considering  an  interaction 
impedance.  For  a  gyrotron  (i.e.,  a  cavity  structure), 


P  .Jili 

Pc  QKl 


(35) 


where  P  is  power  in  the  cavity,  Q  Is  the  loaded  Q,  E  is  the  rf  field  strength,  and  K_  is 
C  L 

related  to  the  interaction  impedance.  One  need  only  to  determine  for  the  particular 
cavity  in  question.  With  the  power,  efficiency,  and  d.c.  voltage  (which  determines  the 
value  for  v„),  the  optimum  operatiag  current  is  obtained  from 


> 


I  -  PCA,VA. 


(36) 


In  summary,  the  following  optimum  operating  point  characteristics  of  a  particular 
design  can  be  obtained  from  the  nonlinear  theory! 

d.c.  magnetic  field 
rf  field  strength 
rf  output  power 
efficiency 
current 
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where  the  input  parameters  are  the  frequency,  the  d.c.  voltage,  y* ,  v„  ,  and  the  parallel 
and  perpendicular  wave  numbers.  The  only  additional  parameters  necessary  are  those  to 
convert  rf  field  strength  to  output  power. 

4.0  COMPARISONS 

4.1  COMPARISON  OP  NONLINEAR  THEORY  AND  COMPUTER  SIMULATIONS 

In  order  to  verify  that  the  nonlinear  theory  is  correct,  a  comparison  is  made 

between  the  optimum  operating  characteristics  for  a  particular  design  as  predicted  by  the 

nonlinear  theory  versus  those  obtained  from  a  computer  simulation.  The  particular  design 

used  for  this  comparison  is  as  follows; 

f  -  60  GHz 
o 

k„  =  1.013  cm-* 
kA  =  11.14  cm-* 
vw  *  0.335  c 
vA  *  0.45  c 
Va  *  80.0  kV 

With  (31)  through  (36),  the  following  are  obtained  from  the  nonlinear  theory; 

B_  -  33.02  kG 
P 

Ef{  »  2.43  X  10 7  V/m 
Pc  -  324  kW 
n  *  58.6796 
I  -  6.9  A 
v*#  -  0.294c 
avA  -  0.034c 
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The  loaded  Q  for  this  particular  design  is  400  and  which  is  determined  from  a  Varian 

o 

computer  code  that  solves  the  eigenvalue  problem  for  the  cavity  is  4.545  X  10  yielding  P 

c 

-  324  kW  for  Erf  -  2.43  X  10?V/m. 

The  results  from  the  computer  simulation  for  the  above  design  are  shown  in 
Figure  4,  (Shively,  et.  al.,  1982).  The  computer  code  that  was  used  is  a  two  dimensional, 
fully  relativistic,  large  signal  code  developed  at  Varian  Microwave  Tube  Division  that 
includes  both  Bff  and  with  no  averaging  over  a  gyroradius  where  the  effects  of  the 
electron  beam  on  cavity  Q  and  the  rf  field  profile  are  not  included.  In  Figure  5,  the  point 
in  parameter  space  obtained  from  the  theory  is  shown  with  the  results  of  the  computer 
simulation.  Consider  the  following:  First,  the  theoretical  value  of  efficiency  at  6.9A  lies 
very  close  to  the  value  obtained  with  the  computer  simulation;  second,  even  though  the 
theoretical  optimum  point  in  parameter  space  is  higher  in  current  than  the  optimum  point 
from  the  computer  simulation,  (l.e.,  6.9A  versus  5.7A),  this  discrepancy  is  most  likely  a 
result  of  incomplete  simulation  results.  The  most  information  obtained  from  Figure  4  is 
that  the  optimum  operating  point  is  where  the  magnetic  field  is  between  22.93  kG  and 
23.05  kG  and  the  current  is  between  5.7A  and  8.8A.  Clearly,  the  optimum  point  is  closest 
to  23.05  kG  and  5. 7 A.  It  is  possible  that  the  optimum  operating  point  is  the  same  as  the 
theoretical  value--23.02  kG  and  6.9A.  In  Figure  5,  the  dashed  line  indicates  a  possible 
simulation  curve  for  a  magnetic  field  of  23.02  kG.  Such  a  curve  indicates  that  the 
nonlinear  theory  accurately  predicts  the  optimum  operating  point  for  a  particular 
gyro-device  design. 

In  addition  to  comparing  the  theoretical  and  computer  results  for  the  optimum 
operating  point,  it  is  useful  to  make  comparisons  at  other  values  of  magnetic  field  and 
current.  Theoretical  predictions  are  obtained  from  (32)  by  substituting  various  values  of 
for  v4j  -  a v^ where  v'>^.  The  theoretical  values  for  v*»  0.325c  and  0.294c  are  23.34 
kG,  5.08A  and  23.18  kG,  6.17SA,  respectively.  These  values  are  compared  to  the 


computer  results  In  Figure  6  demonstrating  close  agreement.  From  the  nonlinear  theory, 
the  line  shown  in  Figure  6  is  generated  with  the  endpoint  being  the  optimum  operating 
point. 

4.2  COMPARISON  OF  THE  NONLINEAR  THEORY  AND  EXPERIMENTAL  DATA 
Having  shown  close  agreement  between  predictions  of  the  nonlinear  theory  and 
results  of  the  computer  simulations,  a  comparison  is  made  between  predictions  of  the 
nonlinear  theory  and  available  experimental  data.  The  experimental  data  (Figure  7)  were 
obtained  from  a  Varian  60  GHz  gyrotron  operating  under  pulsed  conditions.  From  Figure  7 
(Felch,  et.  al.,  1982),  the  following  experimental  values  are  obtained; 


I(AMPS) 

°,E 

f  (GHz) 
o 

nE(%) 

PE(kW) 

B^kG) 

4.0 

23.30 

59.68 

23.5 

75 

24.0 

6.0 

23.17 

59.67 

27.2 

135 

24.1 

8.1 

23.07 

59.67 

32.5 

210 

24.2 

The  values  of  vA  for  the  various  currents  are  obtained  from  the  following; 

bhE  =  bh  <Ey  "  °)  lv±  (37) 

Assuming  that  the  values  of  Bp  given  in  the  preceding  table  satisfy  (32)  and  using  the 
values  of  v,  obtained  from  (37),  ^  and  are  determined.  With  the  cavity's  obtained 
from  a  Varian  computer  code,  the  rf  field  is  converted  to  power.  The  efficiency  is 
obtained  by  substituting  v^  and  vx  into  (33).  The  complete  results  from  the  theory  are 
summarized  below; 


I(AMPS) 

vi 

PT(kW) 

tiT(96) 

4.0 

0.420c 

70.9 

30.75 

6.0 

0.432c 

160.1 

43.60 

8.1 

0.440c 

234.5 

51.49 
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Comparing  the  theoretical  results  with  the  experimental  data,  a  number  of  points  can  be 

made.  First,  considering  the  gun  design,  the  values  of  a(a>vA/vH)  ere  reasonable  (l.e., 

1.87<a<1.96).  Second,  the  power  predicted  by  the  theory  is  well  within  20%  of  the 

experimental  values.  Third,  the  theoretical  efficiency  is  significantly  higher  than  the 

experimental  values,  but  the  computer  simulation  results  are  no  better.  Of  course,  both 

the  nonlinear  theory  and  the  computer  simulation  assume  no  velocity  spread  in  the 

electron  beam.  It  is  likely  that  velocity  spread  effects,  as  well  as  effects  due  to  the  fact 

that  all  the  electrons  are  not  located  at  the  electric  field  maximum  of  the  TB  *  mode 

o2 

account  for  the  lower  efficiencies  seen  in  the  experimental  data.  In  light  of  these 
caveats,  one  can  conclude  that  there  is  reasonable  agreement  between  the  nonlinear 
theory  and  the  experimental  data. 
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5.0  CONCLUSION 

It  has  b«en  shown  that  the  soliton-based  nonlinear  gyro-device  theory  allows  one 
to  predict  the  optimum  operating  point  for  a  particular  design.  The  optimum  current, 
magnetic  field,  efficiency,  and  power  are  obtained.  Comparisons  of  the  theoretical 
predictions  to  computer  simulation  work  done  at  Vartan' s  Microwave  Tube  Division  show 
close  agreement  not  only  for  the  optimum  operating  point,  but  for  other  points  in 
parameter  space  as  well.  Comparisons  of  the  theory  with  Varian  60  GHz  gyrotron  data 
show  reasonable  agreement  when  the  effects  of  velocity  spread  and  electrons  off  the  B 
field  maximum  are  considered. 

Having  demonstrated  that  the  nonlinear  theory  agrees  with  simulation  and 
experiment,  work  can  begin  on  understanding  the  full  effect  of  the  nonlinearity.  Some 
areas  of  immediate  interest  are  harmonic  generation  and  modulational  effects  (i.e., 
amplitude,  frequency,  and  phase). 
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FIGURE  2  HIGHER  ORDER  ELECTRON  MOTION 
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